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O ■ Abstract 

o _ 

We construct majorant functions z) for the Mayer series of pressure satisfying a non- 
linear differential equation of first order which can be solved by the method of characteristics. 
\ The domain \z\ < (er) _1 of convergence of Mayer series is given by the envelop of charac- 

Qh' teristic intersections. For non negative potentials we derive an explicit solution in terms of 

the Lambert VF-function which is related to the exponential generating function T of rooted 
trees as T(x) = — W(— x). For stable potentials the solution is majorized by a non negative 
potential solution. There are many choices in this case and we combine this freedom together 
with a Lagrange multiplier to examine the Yukawa gas in the region of collapse. We give, in 
this paper, a sufficient condition to establish a conjecture of Benfatto, Gallavotti and Nicolo. 



> 



For any f3 £ [4ir, 8tt), the Mayer series with the leading terms of the expansion omitted (how 
many depending on (3) is shown to be convergent provided an improved stability condition 
holds. Numerical calculations presented indicate this condition is satisfied if few particles are 
involved. 



1 Introduction 

One of main investigations in equilibrium statistical mechanics concerning the determination of 
phase diagrams is to define regions of the parameter space in which a phase transition does not 
occur. According to Lee- Yang Theorem |YL| ILYj theses regions are free of zeros of the partition 
function which, for lattice models with attractive pair potentials, are located in the unit circle 
of the complex activity plane. Another way of determining whether the phase diagram is free of 
singularities is given by Mayer series. Although the theory based on these series does not describe 
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the equation of state beyond the condensate point, it is able, however, to deal with continuum as 
well as lattice space for a larger class of admissible pair potentials. 

The domain of convergence of Mayer series for classical gas with stable two-body summable 
potential can be established using various different techniques (see Ruelle jEj, Brydges [B] and 
references therein). There are cases, including attractive potential at short distances and Yukawa 
gas at inverse temperature (3 up to An, which require the use of finitely |GM| or infinitely |Be[ H] 
many Mayer series iterations. The present investigation uses the flow equation of Ursell functions 
introduced by Brydges and Kennedy jBKl IBKlj in order to obtain previously known and some new 
results on the convergence of Mayer series using methods of partial differential equations (PDE). 

The radius of convergence of a Mayer series is often determined by a "non-Physical" singularity 
in the complex activity z-plane. In the present article we shall write a majorant series $(t, z) 
for the pressure p(/3, z) of a gas which satisfies a nonlinear partial differential equation of first 
order. The "non-Physical" singularity can be recognized as the location, in the phase space Q = 
{{t,z) G M 2 : t G [to,tx] , z > 0}, of characteristics intersection (t is eventually related with (3). As 
in Tonks' model of one-dimensional rods with hard core interaction, the singularity of $(t, z) is 
manifested by the W^-Lambert function, an indication that this singularity is of combinatorial 
nature jKj. 

In the application to the Yukawa gas, the variable t G [to, 0] parameterizes the decomposition 
of the Yukawa potential into scales, with t being related with the short distance cutoff, which is 
removed if to — > — oo. There is a conjecture that the radius of convergence of the corresponding 
Mayer series remains finite as to ~~ > — °o f° r P £ [47r, 8n) if one subtracts the leading terms of 
the expansion, how many depending on (3. Stated as an open problem in |Bej and previously 
formulated in |BGN| IGNj . the conjecture will be referred here as BGN's Conjecture. Brydges- 



Kennedy |BKj have proved convergence of Mayer series for 4tt < f3 < 167r/3 with the O (z 2 ) term 
omitted and shown how it could be extended up to (3 < 6tt if an improved stability condition for 
three particles holds. In the present work we establish BGN's conjecture up to (3 < 8tt, assuming a 
mild improvement on the stability condition for the energy of non-neutral n-particle configurations. 
As our calculation indicates, this condition can be easily verified numerically if n is not too large 
since the minimal value attained by the energy maintains a comfortable distance from the lower 
bound needed to establish convergence. 

At inverse temperature 8ir (1 + (2r + , r = 1,2, . . ., neutral multipoles of order < 2r 

become unstable and collapse if the cutoff is removed. Even coefficients of the Mayer series up to 
order 2r diverges as to ~~ > — 00 an d they remain finite otherwise because there occur cancellations 
when unstable multipoles interact with another particle. Our proof of convergence has two distinct 
parts. In the first, the solution & k \t, z) of a nonlinear first order PDE modified by Lagrange 
multipliers, with $( fc )(to,z) = z, is shown to be analytic in a disc whose radius remains positive 
as to — > — oo for any (3 < (3k = 8tc (1 + A; -1 ) . In the second part, $^ fc - ) is shown to majorize 
the Mayer series with even coefficients up to order k subtracted. The two ingredients here are an 
improved stability bound for non-neutral configurations and translational invariance. All estimates 
are done in the infinite volume limit to avoid boundary problems. As the Yukawa potential decays 
exponentially fast this choice, although not relevant from the point of view of thermodynamics, is 
suitable to exhibit the necessary cancellations right at the beginning. 

This paper is organized as follows. Section |2] contains a review of Gaussian integrals, including 
their relation to the partition function of classical gases, and the statements of our results in 
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Theorems 12.21 and 12.31 Section |3] gives an alternative derivation of Ursell integral equation of 
Brydges-Kennedy. The integral equation is the starting point for the Majorant differential equation 
presented in Section This first order nonlinear partial differential equation is solved in Section 
|3] by the method of characteristics for non negative and stable potentials, concluding the proof of 
Theorem 12.21 Section |U] and [7| give an application of the previous section to the ultraviolet problem 
of the two-dimension Yukawa gas. Section El reestablishes previously known results via our method. 
In Sections [3 we prove BGN's Conjecture under the assumption of an improved stability condition 
(Theorem 12. 3j) which has been verified numerically for n = 3,4 and 5 (see Remark l7.3j) . 



2 Gaussian Integrals and Summary of Results 

To fix our notations, we need some facts on Gaussian integrals and their relation to the grand 
partition function (see references |FSt IBr| iBMj . for details). 

The state of a particle £ = (x, a) consists of its position x and some internal degrees of freedom 
a. Examples: 1. For one-component lattice gases, ( = iG Z d ; 2. For a two-component continuous 
Coulomb gas, £ = (x, ±e) where x G M. d and e is the electrical charge; 3. For dipole lattice gas, 
£ = (x, ±ej) where x G Z d and ej is the canonical unit vector. 

The particles are constrained to move in a finite subset X of Z d (or a compact region of M. d ). 
The lattice is preferred for pedagogical reason but we shall also consider particles in a continuum 
space. S is in general a compact space and we set A = X x S. By the triple (Q, E, dp) we mean the 
measure space of a sing le particle with Q = R d x S , E a sigma algebra generated by the cylinder 
sets of M. d x S and dp(£) = d/i(x) dv(a) a product of positive measures (dfi(x) = ^ d 5(x — n) d d x 

is a counting or Lebesgue dfi(x) = d d x measure times a probability measure dv on S). 

The grand canonical partition function is given by 



OO yi p 

n=0 H - Jkn 

where, for each state £ = ...,£„) of n-particles with activity z, 

Mi) = z~ PUn{i) (2-2) 

is the associate Boltzmann factor, 

0»(*)=X>(&6) (2-3) 

i<j 

the interaction energy and (3 is the inverse temperature. The two-body potential $ : A x A — > M. 
is a jointly measurable function and we use, sometimes, the abbreviation := The finite 

volume pressure pa(/3,z) of a classical gas is defined by 

PA(M = ^hS A GM. (2-4) 

Without loss, is defined in the whole space Z d (or M. d ) - the Green's function of some differential 
equation with free (insulating) boundary conditions - and it is assumed to be of the form 

/M(£,0 = aa'v(x-x') (2.5) 
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for some real-valued measurable functions v satisfying: 



1. Regularity at origin \v(0)\ < 2B < oo; 

2. Summability 

dfi(x) \v(x)\ < oo ; (2.6) 



3. Positivity 

2j Zi Zj v (xi - Xj) > , (2.7) 

l<i ,j<n 

for any numbers x±, . . . ,x n E M. d and Zi,...,2„ 6 C, n = 1,2, .... 
Note these assumptions ensure that •& is a stable potential, i.e., choosing z\ — • • • = z n = 1 in 

(EH), 

>-i E ^K0)l>--Bn (2.8) 

l<i<n 

holds with B > independent of n (see Section 3.2 of Ruelle [H], for further comments). 

Since v is a positive definite symmetric function, it defines a Gaussian measure \x v with mean 
dfi v ((f>)(f>(x) = and covariance J dfi v (cf)) 4>(x) cf)(x') = v(x — x'). Introducing the inner product 



{fay) = J dfi(x) (f)(x) r](x) 

and the indicator function 

n 

yn(-) = y(Z;-) = J2 a i s (-- x ^ ( 2 - 9 ) 

3=1 

of a n-state £ = (£i, . . . , £ n ), we have 



^n(0 = y dfrte) :e^: v (2.10) 
where : :„ is the normal ordering with respect to the covariance v: 

: e ia ^ x) : v = exp |y w (°)} ■ (2- 11 ) 

Substituting (pHUj) into flTIJ) yields 

E A (/3, z)= f dp* (0) exp {V o (0)} (2.12) 



where 

V (<f>) = z I dp(0 :e^: v . (2.13) 



A 
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Example 2.1 1. V (<f>) = z ^ : e i4>[x) : v for one-component lattice gas; 

2. Vq{4>) = z d d x : cos0(x) : v for standard Coulomb or Yukawa gas; 

Jx 

3. V o (0) = z 2^ / d d a : cos (a ■ V0(x)) : v for dipole lattice gas. 

In Example 2 we have respectively, v = (— A) -1 and v = (—A + 1)~ , the fundamental solution 
of Laplace and Yukawa equations. The potential of two dipoles, of moment a and a', in Example 
3. is given by v = (er • V)(cr' ■ V')(— A) -1 , second difference of the fundamental solution of discrete 
Laplace equation. Vo in Examples 2 and 3 contain cosine function instead exponential. This 
form expresses the pseudo-charge symmetry required for the existence of thermodynamic limit 
lim A ^ R d( Z d) pa —p of respective systems (see, e.g., jS]). 

The potentials in Examples 2 and 3 do not fulfill all requirements: Coulomb and dipole potentials 
do not satisfy summability whereas the Yukawa potential fails to be regular at origin. To circumvent 
these problems, a one parameter smooth decomposition of covariance v(t,x) may be introduced 

so that v(to,x) = and / \\v(s, -)|| ds < oo holds for any finite interval [to,^ij- This allows us 

to approach the ultraviolet problem of Yukawa gas by letting i ~~ > — oo maintaining t\ fixed (e.g. 
t\ = 0) in Example 2, and the infrared problem of dipole gas in Example 3 as the limit t\ — ► oo 
with to = fixed. The application in the present work deals only with the former problem. 

In all examples Vo(<j>) is a bounded perturbation to the Gaussian measure \i v satisfying |Vo(0)| < 
z\A\. Brydges and Kennedy [BK j have used the so called sine-Gordon representation (j2.12|) to 
derive a PDE in the context of Mayer expansion. Although based in their work, our PDE approach 
in the present article is remarkably different. 

It is convenient to write a flow 

V A (t, z- 4>) = In (/^ } * exp { W)}) (2.14) 

starting from Va (0, z; 0) = Vq(</>) where \i v * F denotes convolution of F by the Gaussian measure 
fi v which, by the Wick theorem, is given for the one-component lattice case by 

/<WCW + = eJi £ v{ t, x - y) ^A FW . (2.15) 

I x,y£Z d I 

To derive the standard Mayer expansion, the covariance v(t) is chosen to be any ( tf 1 parameterization 
such that f (to) = and v(ti) = v. The parameterization involving multiscales, v(t) is required to 
be of the form (|2.5jl satisfying properties 1 — 3 with B — B(t) < oo for all t G [to, ti]. v(t) is said to 
be admissible if fulfills all these properties. Note that Vo(4>) also depends on t through the normal 
ordering : :„(*). 

Equation (J2.14|) can be written as 



e 



E A ((3,ze^) (2.16) 
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where Ha(/3,io) generalizes the grand partition function by including a complex valued activity 
w = w (0 which depends on the particle state: 

oo „ 

E A ((3, ™) = ]T - / d n p{Z) w n {£) MZ) (2-17) 



with w n (£) = Y\_ w (£j) anc ^ Q^.IJI is recovered if = z. 
We also have 



P \ PA ((3, z)\ = |- |VX (ti, z, 0)| < sup |Va (*!, z, 0)| (2.18) 

|A| |A| 

and, as t goes to to (the infinite temperature limit), 

lim V A (t, z, 0) = Vo(0) = z|A| • 

t— >to 

The flow {Va(£, z, 0)/ |A| : < t < 1} provides us an upper bound for the pressure of gases described 
above, starting from the pressure of an ideal gas. 

The results presented in this article can be summarized as follows. 
Theorem 2.2 Let $ = Q(t,z) be the classical solution of 

<S> t = h{z® z f + B{z® z -z) (2.19) 

defined in to < t < t\, z > with Q(to,z) = z where T = T(t) = \\v(t)\\, B = B(t) is defined by 
with v replaced by v(t). Then 

holds uniformly in A. The solution $ has a power series 

oo 

®(t,z)=Z + J2 A n(t) ^ 
n=2 

with A n (t) > 0, convergent in the open disc D = [\z\ < (er) -1 } where 

r(t ,t) = inf QT r(s)exp|2^ B(s') ds'^ds, exp B{s') ds'J J Y{s)ds 



and is majorized by 



<$>(t,z) < — (w(-zt) + Uv 2 {-zt) ) (2.20) 



where W(x) is the Lambert W -function defined as the inverse of f(W) = We . 
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Theorem 2.3 Given an odd number fceN, suppose an improved stability condition 

Un{t;£)>-{n-5 n )B{t) (2.21) 

n 

holds for all n-particle states £ = (cr, x) such that dj ^ with 



i=i 



, 1 

n 



and n < k. Then, the pressure p(j3,z) of Yukawa gas, with even Mayer coefficients up to order k 
subtracted, exists and has a convergent power series provided 

a a All 
p < (Jk and \z\ < 



where flk = 8tt /(l + k x ) is the k-th threshold. 



Theorems 12 . 21 and 12 . 31 will be proven in Sections and respectively, by majorizing the solution 
of (|2.19|) . possibly modified by a Lagrange multiplier L, by an explicitly solution of (|2.19|) with B = 
and T replaced by T/f for a conveniently chosen function / . 



3 Ursell Flow Equation 

The power series of the pressure (12. 4|) in the activity z is called Mayer series: 

oo 

l3p A (J3,z) = ^2b Atn z n . (3.1) 

n=l 

The Mayer coefficient b A>n has an expression analogous to (|2.1|) . 

|A|&A,n = ^ / d n p{d)^ n {t) (3.2) 

n\ J A n 

where, for each state £ = (£i, . . . , £ n ) e (R x S) n of n-particles, ^n(^) i s ^ e Ursell functions of 
order n. In terms of Mayer graphs, 

= £ II (c-^-l) (3-3) 

G (iJ)6G 

where G runs over all connected graphs on /„ = {1, . . . ,n} (G is connected if for any two vertex 
k, I 6 I n there is a path of edges E I n x I n in G connecting k to I). 

Equation ()3.3j) is not useful to investigate whether (|3.1j) is a convergent series since E G 1| > 
c (n!) 1+ holds for some c, 5 > 1 . We shall instead derive a flow equation for ^ na ^ incorporate 

the cancellations required to reduce the cardinality of the sum in ()3.3|) to 0(n\). 

1 The number of Mayer graphs with n vertices is 2™( n_1 )/ 2 and the number of connected Mayer graphs is of the 

n 2 3n 

same order. For this, note that a tree, the minimal connected graph, has exactly n— 1 bonds, remaining — — + 1 

bonds to be freely chosen to be connected. 
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Given £ = (£j) 3 ->i G (R x S)°°, let P n £ = (£i,...,£ n ) denote the projection of £ into the 
subspace (R x S 1 )™ generated by its n first components. If / C N is a finite set / = {ix, . . . , i n }, we 
set Pj£ = (£ il5 . . . ,£i n ). To each sequence £ we associate a complex Boltzmann function 



u(t, £) = exp < - (7i(7jv(t, Xi - xj) + z <?j4>{xj) 



i<3 

and define 

u n (t, £) = u(*, P n £) = e -^n«)+<(^) , (3.4) 

for n > 2 where U n (£) = U (P n £) is the interaction energy ()2.3j) and f] n (-) = r](P n £, •) is the indicator 
function (|2.9jl . Note that w n : R + x (R x S 1 )™ — > C is symmetric with respect to the permutation 
7r of the index set {1, . . . , n}: 

U(t, P n £) = u(t, M n P n £) 

where M n : . . . , £ n ) — > . . . , £ nn ) is the permutation matrix corresponding to 7r. Note also 
that u n (t, £) can be written as 

«n(*, = MS) ^' Vn) = ^(t)* : e i( ^ :«( t ), (3.5) 
where ?/v is the Boltzmann function (|2.2jl . with u = t> (t) being an admissible covariance. 

Now, let u(t) = (it n (t, £)) n>0 denote the sequence of these functions with uo(t) = 1 and Ui(t, £) = 
e tt7l ^ Xl \ At t = to, v(to) = (absence of interactions) and it(to) = u ^ is the sequence with 



u 



(0) (£) = e^vn) f or n > i an d v$\g) = I. Using (jS3J), (l2~TT|) and (j2HSJ), a simple computation 



yields 

Proposition 3.1 The sequence it(t) = (u n (t,£)) n>0 is the formal solution of the heat equation 

u t = -A^u (3.6) 
with initial condition it (to) — it 1 - ^, where the action of Ad on u is a sequence A^u with 

{A D u) n = V D{x hXj ) ®\l n ,. (3.7) 



Remark 3.2 1. The derivatives in \3. T\) are partial derivatives with respect to <fi at the value x, 
even for continuous space. This has to be contrasted with the functional derivative in h2.1ty 
whose correspondence to partial derivative can only be made for lattice space. 

2. To emphasize the order of operation, we write Ap = V • D V for the differential operator \3. 7| ) 
if D is a off-diagonal symmetric matrix. 

Let 6 denote the vector space over C of sequences t/> = (V ; n(€)) n >o °f bounded Lebesgue 
measurable functions ip n {£.) = ^{Pn £) on (R x S) n which are symmetric with respect to the per- 
mutation 7r of the index set {1, . . . ,n}, with ^>o(£) — 4> (P$£) £ C. We define a product operation 
o : © x © — > & for any two sequences of symmetric functions ip 2 as a sequence 

1> 1 Oll> 2 = ll>=(l(> n (£))n ( 3 - 8 ) 
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whose components are given by 

IC{l,...,n} 

where I c = {1, . . . , n} \J. Note that this operation is well defined and preserves the symmetry. & 
endowed with the operation o form a commutative algebra with the identity 1 = ) n>0 given 

by 

1 if 1 = 
if 1^0 ' 



Let ©o and &i denote the subspace and the affine space of sequences if) = C0n(O)n>o such 
that Vo(£) — and ^o(€) = 1> respectively. We define the algebraic exponential of sequences 

4> = (Mt))n>0 

Sxp : &o — > &i 

as the following formal series 

Sxp (0) = 1 + + ^0 O (/) + ■■■ + —(() O ■ • • O + • 

2 n! v 



n— times 



and the algebraic logarithmic of sequences tf) = {i> n {£)) n >o 

£>n : &x — ► 6 

as (ip e ©o) 

1 f-D n_1 
J&fn (1 + V>) = t/> - o if) -\ h - — V ° • • • o i/H . 

n— times 

Note that the following 

Sfn(£xp(<f>)) = 4> and #a;p (jgfn (1 + V)) = 1 + 1> 

holds for any <f>, if} E & , so <fa;p and J5fn are algebraic inverse function of each other. We thus 
have (for details, see Ruelle ITT] ) 

Proposition 3.3 The sequence i/> = (V ; n) n>0 °/ Boltzmann functions given by n > 2 mi/i 

= V'i(^) = 1> an algebraic exponential of the sequence if) c = (V , n) n >o °/ Urs ell functions given 
by (O wit/j = and ^ c (£) = 1; 

i/; = <f (V> c ) ■ (3.9) 
Proof. Let (•, •) : 6 x 6 — ► C be the inner product 

00 n p 

M) = Y. Z - <t>n($) Mt) (3-10) 

n=0 ' J 
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and let x = (Xn(0)n>o ^ e ^ ne sec i uence °f characteristic function of the set A": 

n 

x»(€) = IIxa(^) 

with 

m = J 1 if * G A 

X ™ \ otherwise ' 

Using 

(x,0ot/>> = (x,<f>) (x,i>) 

((f> — > (x, <f>) = (x, defines an homeomorphism of & into the algebra of formal power series) 

together with flH]|, lETTfll) and (JXTJ), we have 

SaGM = (x,V>(*) = (X^xp(^ c ))(z) = ex P {( X ,^ c )(z)} 
which, in view of (|3.1|) - (|3.3J) . proves the proposition. 

□ 

Replacing if) in (|3.9|) by the complex Boltzmann sequence u(t) = (u n (t, £))n>o> l ea ds to the 
complex Ursell sequence y(t) = (yn(t,£)) n>0 given by 

y(*)=J&?n(u(f)) . (3.11) 

Note that, in view of (EHTl) . (l3~oT) . (l3~THl and the fact 

(u o u ) n (t, $ = J2 w> ^ p i c & e ^' Vn) = u> ° vo„ & > 

JC{l,...,n} 

we have y n (t,£) = Vn(*>£) e ■ We shall use ()3.11|) together with ()3.6|) to derive a partial 
differential equation for the complex Ursell functions y n (t,£). A differential equation for the Ursell 
function if>^(t,^), given by (|3.3j) with f3i9(£i,£j) replaced by ai(jjv(t obtained setting; 

= at the end of computations. 

The derivative of (|3.1H1 together with Remark 13.21 2. yields 

1 

y = — o U t 

= -Sxp(-y)o(V ■vVSxp(y)) 

= -A*y+-Vy«i)Vy (3.12) 



where 



(Vy.DVy) n = ^ D(x t ,x j )( 

1 i-l-i V 



dy dy 

o 



l<i ,j<n,i^=j 



d<p(xi) d(f)(xj)J n 



£ ^y{Pii)-DVy(P^) 

/C{l,...,n} 
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from which equation (6) of |BKlj is obtained: 



J^(£) = -Un® wt) - \ E u{p&Pi4)r{Pi<i)r{Pi4) (3.13) 

JC{l,...,n} 

satisfying initial condition t/>q = (^o, n (^)) n > a ^ with 

*(^ = {i othe^se ' < 314 » 

For brevity, we have omitted the time dependence in 1)3.13)1 . Here, U n (t,£) = U(t,P n £) is the 
energy of state P ra £ and 

U(t,P I &Pj£) = £ rr ; rr / r ( /.,. . (3.15) 

The minus sign in the r.h.s. of (j3.13|) comes from the derivatives of at = 0. 

Using the variation of constants formula, ()3.13j) is equivalent to the following integral equation 

If* - / dTU n (r,£) . 

r n (t,t) = -~ dse J* u^p&p^r&Pitw&Pi't), (3.16) 

*° /C{l,...,n} 

for n > 1 with %{t, £) = and ^0, 1) = 1 for all t G [t , 
The second Ursell function ip% satisfies 

rt 

>•* -0-10-2 / dr v (r, X\ — x 2 ) 



f 

f 1 —<Tl<T2 I 

V'lfo&jfa) = ~ dsa 1 a 2 v{s,x 1 -x 2 )e Js 

Jtn 



't 

-01CT2 / dr-?) (r, xi — x 2 ) 
= e J *o -1 (3.17) 

and agrees with (J3.3)) if ^(^1,^2) — ( v {ti, %i — ^2) — ^(^oj^i — ^2))- One can verify that 
there is a unique solution to the integral equation (J3.16)) formally given by the Ursell functions 
(|3.3j) . These follows from ()3.13|) by induction in n. Equation ()3.16|) . derived for the first time by 
Brydges-Kennedy |BKj . is the starting point for the procedure of the next section. 

4 Majorant Differential Equation 

An one-parameter family g(t, z) of holomorphic function on the open disc 

D R = {z EC:\z\ < R} 

has a convergent power series 

00 

g(t,z) = J2C J (t)z> 
3=0 

which we denote by C • z with C = {Cj)- >Q and z = (V) J>0 . 
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Definition 4.1 g(t, z) = C(t) ■ z is said to be a uniform majorant of a function f(t, z) = c(t) ■ z 
in Q) = [a, b] x D R if 

\ Cj (t)\ < Cj(t) 

holds for every j G N and t G [a, b] . We write f < g for the majorant relation. 

Note that a majorant g is a real analytic function and the majorant relation is preserved by 
differentiation in z and integration with respect to t. 

Proof of Theorem \2. 6 A Part I. Our aim is to find an equation for a uniform majorant Q(t,z) = 
A{t) ■ z of the pressure 1)3.1)) . Let a(t) = (a n (t)) n>1 be given by 

a n (t) = sup / dp(£ 2 ) • • / dp(£ n ) (4.1) 

and define 



r(t):=||t>(t)|| , 7M):= y B(r)dr, (4.2) 

where I? and ||-|| are given by ()2.8)) and (|2.6j) . It follows from 1)3.16)1 that 



(*) < g y ^e^r( S ) ^ k(n-k) a k (s) a n _ k (s) 



(4.3) 



holds for n > 1 with oi = 1. We thus have 



|&A,n| < -lOn(*l) < Mh) (4.4) 

uniformly in A where b^ n is the Mayer coefficients ()3.2)) and A{t) = {A n (t)) n>1 is a sequence of 
continuous functions so that (n\A n (t)) n>1 satisfies ()4.3j) as an equality: 

A n (t) = - ds e n ^T(s) k Ms) {n - k) A n _ k (s) (4.5) 
1 Jt ° k=i 

for n > 1 with A\ — 1. 
Writing 

oo 

$(M) :=^A n (t)z™ (4.6) 

n=l 

equations (jSHJ), (IP)) and (1431) yields 

/3j>aOM <*{h,z) (4.7) 
uniformly in A where $ (ti, 2) is the solution of 



If 1 ( d$ , 

$(t,z)=z + - y rfsr(s) [z-fo(s,. 



ze 1 ^) ) (4.8) 
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at time t = t\. Here, the following identities have been used 

oo n— 1 oo n— 1 



J2z n J2kA k (n- k)A n _ k = ^kA k z k (n-k)A n ^ k z r 



n=2 k = l n-l=l k=l 

2 



Y^kA k z k (n-k)A n ^ k z n - k = (z—j 

k=l n-l=fc - ' ' 



The partial derivative of $ with respect to t and z are denoted by $ t and $ z . From ()4.6|) and 
(Oil , we have 

*e 7 )) t = |(^ + 2Ve 2 ^ 2 ( S ) + 3VeH(3) + -) 
= Bse* ze^) , 

with 7 = 7(5, i), and 

z<t> z = z + [ dsT{s)z<!> z (s,ze^ s ' t) ) ze^ {z<5> z ) z (s, ze^) . 
Differentiating ()4.8|) with respect to t, 



thus yields 

i] 
2 

with initial condition $>(t Q ,z) = z. This establishes the majorant relation of Theorem 12.21 and 
equation ([2.19)1 . 

□ 

It will be convenient to consider the density 

p A (P,z)=z^(P,z) (4.9) 

whose corresponding Mayer series has the same radius of convergence of the Mayer series of p\. 
Writing 9(t, z) = $ 2 (t, z), we have 

-Pa{P,z)<Q(P,z) (4.10) 

z 

with satisfying the following quasi-linear first order PDE 

& t = r (z& 2 + z 2 ee z ) + b (z& z + 0-1) (4.11) 

with initial condition Q(t ,z) = 1. 
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In view of (14. 6|) . ()4.9|) and ()4.1U|) . we seek for solutions of (|4.11|) of the form 

oo 

Q(t,z) = l + Y,C n (t)z n - 1 (4.12) 

n=2 

with C n = nA n (t) > and this form is preserved by the equation. 

Remark 4.2 The convergence of Mayer series for an integrable stable potential (v(t) so that 
ti 

\\v(t)\\dt < oo and v(ti) = fid) can be established as in the case of non negative potential 



to 



(B = 0). Writing w = e l{t °^z and A n = e- n ^ to ^A n , we have 

oo 

*{t,z)=Y,Mt) w n :=$(t,w) 



n=l 



with $ satisfying <3> t = (r/2) \ w$ w j . The Mayer series of these systems converges for all ((3, z) 
so that the classical solution $(£, e 7 *-* '* 1 ^!) exists up to t = t\. 

5 The Method of Characteristics 

In this section the convergence of Mayer series will be addressed from the point of view of (j4.11|) . 
We seek for a subset of {{t, z) : t G [to, ti) ,z > 0} in which there exists a unique classical solution 
of the initial value problem. We shall distinguish two cases, depending on whether the stability 
function B(t) is or is not identically zero. 

5.1 Non Negative Potentials 

Let equation (}4~TTj) with B(s) = 0, 

9 t = T (zQ 2 + z 2 QQ z ) (5.1) 

with Q(t Q ,z) = 1, be addressed by the method of characteristics. This equation describes, by 
(I4.1(J|) . a density majorant of a gas with non negative (repulsive) potential v > 0. 

It is convenient change the t variable by r = / T(s) ds and write U(r) = Q(t(r), z(r)). 

Jto 

Equation 

U' = t'e t + z'B z 

together with ()5.1|) yield a pair of ordinary differential equations (ODE) 

U' = zU 2 (5.2) 
z' = -z 2 U (5.3) 
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satisfying the initial conditions U(0) = 1 and z(0) = z . Dividing ()5.2|) by ()5.3|) gives 

dU U 

dz z 

which can be solved function of z 

U(z) = - . (5.4) 

z 

Inserting (J5.4)) into (|5.3|) . leads to the following characteristic equation z' = —ZqZ whose solution is 

z(t) = z exp{-z r} . (5.5) 

Now, 0(t,z) is given by ()5.4|) with z = zq(t, z) implicitly given by (j5.5|) . If W{x) denotes the 
Lambert H^-function jKj, defined as the inverse of f(W) = We w (/ (W(x)) = x and W (f(w)) = 
w), we have 

z (t,z) = — W(-zt) 

T 

and the solution of the initial value problem (j5.1|) is thus given by 

Q(t,z) = —W(-zt) 

ZT 



provided 



ezr = ez T(s) ds < 1 . (5.6) 



to 



Since Q(W(x)) ^ for x < — 1/e in the principal branch of the W -function (see e.g. [K]), this 
condition has to be imposed to ensure uniqueness. W has a branching point singularity at — 1/e + zO 
with W(-l/e) = &(W(-l/e)) = -1. At the singularity z* = z*(t) = (zr)" 1 , Q(t,z*) = e and the 
series (j4.12j) diverges. 

Remark 5.1 1. For non negative potentials we can actually find a pressure majorant Q(/3,z) 
explicitly. Differentiating x = W(x)exp (W(x)) and solving for W gives, after some manip- 
ulations, 



and this, together with H4.y\) and W(0) = 0, yields 



2N ' 



z) = / ^^dx = [W (-ZT) + -W 2 (-ZT 

T Jo X T V 2 

2. The interpretation of \5. 6]) is as follows. Let ^(r, Zq) = (lnr — ln^o) / (r — zo) be the time 
necessary for two characteristics starting from r and zq to meet at a point z: 

z = re~ rSr = zoe~ zoSr . 

^(r,zo) is a monotone decreasing function of r with lim r y Zo £?{r, zq) = 1/zq. Note that, 
Q(^, z) assumes two different values U(z) = r/z and Uq(z) = zq/z for r < zq which tends to 
e in the limit r /* Zq. The minimum intercepting time of two characteristics with r < zq is 
attained when both coincide and this occurs exactly at the branching singularity —zS? = — 1/e 
ofW(x) where uniqueness of \5.1\) breaks down. 
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3. From condition \5. 6}) and J^.7| ) ; the Mayer series \3. 1)) converges if 



\z\ < (e jf r(s) ds^ . (5.7) 

rh 

Choosing v(t) — (t — t )j3$/(ti — t ), we have / T(s) ds = (3 and /or stable potentials, 

___ Jto 
according to Remark \4-fy the Mayer series converges for all (f3,z) so that 

e\z\exp^£ B(s) dsj J V{s) ds <1 (5.8) 

which gives the same condition eft \z\ e^ B ||$|| < 1 obtained by Brydges and Federbush (see eg. 
(11) in JEE})- 

5.2 Stable Potentials 

We now turn to equation (j4.11|) . Instead of solving it explicitly, we shall obtain an upper solution 
0(i, z) whose existence implies convergence of density majorant series. Q(t,z) satisfies equation 
(|5.1j) for non-negative potential with T substituted by T/f for some properly chosen functions /. 
Our procedure allows us to obtain various convergence conditions depending on the choice of /. 

Proof of Theorem \2.°A Conclusion. Let us change the variable z by w = f{t)z, for a function 
f{t) > 1 such that /(0) = 1, and define 0(t, z) = f(t)^f (t, f{t)z). Using the chain rule, we have 

( f f 

zQ z = fw^ w 
which, together with (|4.11jl . yields 

% - \B - 4 + Twv) w^ w = (b~^]^/ + Tw^ 2 - 4 • 

Choosing 

t 



f J V fj f 



fit) = exp <j / B(s) ds 

yields 



t 



- Tw 2 ^ w = Tw^ 2 ~j< Tw ^ 2 ( 5 - 9 ) 



by positivity of / and B. Note that \1/ satisfies (|5.1|) as an inequality. 



As before, we define r = / Vis) ds and write 



to 



U(t) = ^{t{r),w{r)). (5.10) 
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From U' = ^ t t' + ty w w'axid (|5~§j) . we have 

[/' < wU 2 (5.11) 

w' = - W 2 U (5.12) 

with initial conditions U(0) = 1 and u>(0) = wo- Equations (|5.11jl and (|5.12J) can be written in 
terms of the variables V(r) = wU and 10 as 

V < o 

-y (5.i3) 



with V(0) = lOo, which implies 

V(r) < w (5.14) 

and 

ty(r) > w exp{— w r} . (5.15) 

Here, we need an inequality for wq = wo(t, z) so that ()5.15|) holds. Writing ()5.15|) as 

-w t e~ w ° T > -wr , 

we have 

w (t,w) < — W(-wt) , (5.16) 

T 

where W is the Lambert jy-function, provided 



wt = zexp i^J B(s') (is' | J T(s) ds < 1/e . 



Plugging (|5.16|) into (J5.14|) . together with (|5.1U|) . leads to 



9(t ; z ) = YaDSI < 0( t; z ) } (5.17) 



where 



e{t,z) = w(-zexp{[ B(s) dsX [ T{s) ds 

z [ T(s) ds V Ut ° J Jt ° 

recovering the statements of Remark 14.21 

We now go beyond inequality (|5.15j) . For this we modify slightly the previous scaling by defining 
0(t, z) = (t, fi(t)z). By the chain rule, is the solution of 

*t - ^w 2 W w = (b-§-)w^ w + B^-B + ^rf 
J i V hj Ji 

satisfying ^^(t ,w) = 1. Choosing 

hit) = exp 1 2 J B{s) ds\ (5.18) 
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to 



and using 

oo oo 

n=2 n=l 

oo 

= -^2(n-2)C^w n ~ l < 

n=2 

in view of ()4.12|) and w = zfi(t) > 0, we have 

r r r 

^ t - — w 2 tyty w = B (-wV w + V - 1) + — w^ 2 < —w^ 2 . (5.19) 
Ji Ji Ji 

We define r = / T(s)/fi(s) ds and repeat the steps from (|5.1()jl to ()5.17|) . Provided 

•/to 

WT = ZTl < 1/e (5.20) 

holds with 7i given by 

rt 2 / 5(s) (is 

n (t) = I e Js Y{s) ds , (5.21) 

•/to 

the solution of the initial value problem (|4.11jl satisfies 

Q(t,z) < —W{-t x z) . (5.22) 

TiZ 

It follows from (J5.17j) . (|5.22j) and Remark 15.11 1. that (|2.2(Jj) holds and converges either under 
the condition (|5.8j) or under (|5.2(J|) . This concludes the proof of Theorem 12.21 

□ 

Remark 5.2 1. Working directly with the integral equation \3.1b\) . Brydges-Kennedy JBKU 
have established recursively 

a n (t) < Ti n -V- 2 

with Ti and a n given by \5.21)) and Hence, the density majorant of stable potentials 

satisfies 

a n{t) n _! 

As 



71=1 



00 < \n-l 



71=1 

and i/izs, fry tae Taylor series of the Lambert W -function U^ . is exactly \5.2ty) . 

2. Criterion \5.20\) for convergence may be more advantageous than h5.£fy when v(t) is a multi- 
scale decomposition of a potential. As observed by Brydges-Kennedy (see Remarks at the end 
of IBKljp . a poor stability bound (large B{t)) at short distances (t <C 0) may be overcome by 
the smallness of the interaction T = \\v(t)\\ at these scales. 
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6 Application to Yukawa Gas 



The majorant method is applied to investigate whether the Mayer series for the Yukawa gas is 
convergent in the region of collapse. The first and third subsections redo the proofs in [B | IBK] and 
BK for (3 < Ait and Ait < (3 < 167r/3, respectively. After introducing the model we illustrate how 
useful is condition ()5.20|) for (3 < Air. The procedure of Subsection 15.21 is then modified in order to 
overcome the first threshold at (3% = An. In the present section we intend to show the capabilities 
of our majorant method to extend results in jBKj . The proof of Theorem 12.31 will be deferred to 
Section 

6.1 Convergence of the Mayer Series for (3 < 47r 

Particles in this system have assigned two charges o G { — 1,1} and they interact via the Yukawa 
potential 

0(&,6) = °i°2 (-A + If 1 (xi - x 2 ) 

= / e d 2 k (6.1) 

where —A is the Laplace operator in M 2 . Potential (|6.1j) violates condition 1. of (|2.5jl . since 

|^,0I = - Jim / -^-dr) = — lim In (R 2 + l) 

Air H^oo J Q t] z + 1 Air -R-^oo ' 

diverges logarithmically. We replace ^(^1,^2) by o\o%v{$,x\ — x%) given by 

v(t, x)=/3 (-A + K^y 1 (x) = 7^K ( y/tfjj (6.2) 



where K is the modified Bessel function of second kind (see e.g. |GJj p. 126) and k is a monotone 
decreasing scale function such that 

k(— 00) = 00 and k(0) = 1 , 

and decompose it into scales 

v(t, x) = v(s, x) ds 
J —00 

where 



K(V<sj\x\) ■ (6-3) 



(3 k'(s) , 

V(8,X) = == X 

By definitions ()2.8)1 and (|4.2jl . we have 

Proposition 6.1 For eac/i —00 < s < ; ii(s, x — y) is a positive kernel satisfying properties 1 — 3 
of / OOP 

B( a ) = i|«( fl ,0)| = ^-(-lB «(*))' (6.4) 

2 57T 

and 

r( a ) = 11^,011 = ^^^. (6.5) 
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Proof. Properties 1 — 3 of ()2.5|) follow immediately from the expressions ()6.2|) and (j6.1|) . provided 
(16. 4j) and (|6.5|) are established. Applying Cauchy formula for the integral (|6.1|) (with 1 replaced by 
/t(t)) in the x — y direction, gives 

P 



v(t, x - y) = —g ( V K(t) \x - y\ 



with 

By the chain rule 

(g o w)' 



g(w) :-- 



oo g-ui\/fc 2 +l 



dk 



-W 



—wV k 2 +l 



W 



dk = 



w 



w 



-Xd( = h(w 



w 



(6.6) 



(6.7) 



and equation ()6.3j) can thus be written as 

v(s,x - y) 



-8k'{s) 



4-7T KlS) 



s) \x - y 



(6. 



which gives ()6.4|1 by taking \x — y\ = 0. Note that —k'(s) > 0, by assumption. Integrating ()6.6|) 
times w, yields 

1 



g(w)wdw 



o (Jfe 2 + 1 



.3/2 



rffc = 1 



(6.9) 



which implies 



\v(s,-)\\ 



8 d f°° ( r-- \ 
/ Q \\/K[s)r ) znrdr 







Ms 



and concludes the proof. 



□ 



The majorant density Q(t, z) of the Yukawa gas satisfies equation (j4.11|) with B and V given by 
f!6.4|) and ()6.5|) . respectively, and initial data Q(to, z) = 1 where t plays the role of a short distance 
cutoff, which is removed if to ~~ > — oo. 

By (j4.10|) and Remark 15.2121 the Mayer series converges if ()5.20|) holds. We have 

lim n (0) = -8 lim f f ds = lim {l - k^ 1 ^)} = 

provided 8 < 4tt, independently of the scale function k. The Mayer series of Yukawa gas converges 
if 

47T - 8 

8 < 47T and 11 ' 



z < 



Although it has been previously established in [Bl IBKj , our method is independent of the scale 
function k and improves their radius of convergence. 
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6.2 The 2— Regularized Majorant Solution 



The Mayer series for the Yukawa gas between the first and second thresholds, Air < (3 < Qn, 
diverges because the norm a 2 (t) of second Ursell function ^(i, diverges when the ultraviolet 
cutoff is removed (see e.g. |BGN[ IBKj and references therein). We are going to show that the 
regularized density majorant Q^(t, z), defined by Q(t,z) with the singular part of a 2 (t) removed 
by a Lagrange multiplier, exists and implies convergence of the Mayer series for (3 < 167r/3. For 
167r/3 < (3 < 6n, the upper bound A 3 for the norm of ip%(t, £) also diverges and equation ()4.3|) with 
n = 3 has to be modified in order the majorant relation (j4.10|) holds with O replaced by 0( 2 ). The 
modified equation requires an improved stability condition which prevents (3 2 = 167r/3, and more 
generally (3 = (3 2r , r = 1, 2, . . ., where 

(3 n = 87r(l + -) = 8tt^— , (6.10) 
\ n J n + 1 

to be considered a threshold of the Yukawa model (see Remarks on pg. 47 of BKJ). 

The second Ursell function ip 2 (t,£) * s explicitly given by (|3.17)1 . Together with (|fj.2jl . we have 

V#,£i,6) = exp{^|p (^(v^Wr) -ir (v^)r))}-l (6.11) 

where r = \x\ — x 2 \- 

To isolate the singularity of ()6.11|) we replace a 2 (t) = supg lgA / dp{^ 2 ) \<P 2 {t, d> C2) | by 



x l J ^ _r , ,i 



o-i,o- 2 e{-l,l}: 



which distinguish whether the charge product <J\a 2 is positive or negative {p\Oi = — 1 <^ o\ + a 2 = 
0). Using the fact that K (x) is a monotone decreasing function together with 1 — e~ y < y for 
y > 0, we have 

4(t) < (3 J r(iT ( v ^(t)r)rfr-^o(v / ^)r))rfr 

< P 



.«(*) K (^o) 

For o"icr 2 = —1, the first line of ([3.17)1 . v(s,r) > and (|6.8|) . imply 



a~ (t) = S- r dwwhiw) f -!^$exp(^- [ -^-h ( Jk(t)/k(s)w) dr\ ds (6.12) 
2 Jo J tQ k 2 (s) [4:1tJ s k[t) V J J 



- iK-l3J a 1 ' I k(«) k(«„) V <c(t) y 



/3/4ti^ 



/ 1 1 U{t Q )\ 

4n-(3\ K (t) n(t ) V / 



/3/4tt n 
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Note that h(w) is a monotone decreasing function with h(0) = 1 and 

/*oo /*oo 

dw w h(w) = / dw w 2 g'(w) = 2 dw w g(w) 



(6.13) 



by (jfi.7j) . integration by parts and (|fj.9jl . As a consequence, both a 2 and a^" converges as to — > oo 
if f3 < 47r (a J remains bounded for any /3 > 0). If (3 > 4ir, since the main contribution of ()6.12|) 
comes from a neighborhood of w = 0, there exist £ > and a constant C = C(e) > such that 

— e > 1 and 

47T 

ft(t ) \ K (t) ) K>(t) 
diverges as t — > — oo. 

Now we shall write the majorant equation (|4.11j) with the singularity of A 2 removed by a 
Lagrange multiplier. This can be done in several ways. We may define 6^ such that G?(t, 0) = 
for all t > 0, as the solution of 

@ t = r ( z e 2 + z 2 QQ z ) + B {zQ z + - 1) - z (t + 2Be z {t, 0)) 

with 6(0, z) = 1. Note that d9*(t,0)/dt = 9g(t,0) = holds for all t > 0. Another way adopt 
here is as follows. The linear part of equation (j4.11j) is responsible for the singularity of A 2 if (3 > Air 
and we may choose, instead, a Lagrange multiplier proportional to B 

Q t = T (zQ 2 + z 2 QQ z ) + B {zQ z + - 1) - ^zBQ z {t, 0) (6.14) 

(with a quarter of intensity needed to make A 2 vanishes). 

Let us denote the solution of (RHl) by 9 (2) . As in Subsection EH if 6 (2) (t,^) = * (2) (t,w), 
where w = f2(t)z and 

'3 rt 

to 



/ 2 (t)=expj- / B(s)ds} , ((3.15) 

oo 

using (with ^ 2 ) = 1 + ^2c^w n ~ 1 , C {2) > 0) 

n=2 

1 1 I v v I 

--w^ + ^-l- l -w(¥ 2 \t,0)) = -\^{n-l)C^w^ + ^C^-l-\c 2 w 

n=2 n=l 



1 oo 



n - .'5 



it follows that \E^ 2 ) is the solution of 



% ~ yW 2 W w = (b-^w* w + bU-1- \c 2 w^j + yw& 

= B (-~w^w + # - 1 - + J 2 W ^ 2 

< — w^ 2 

J2 
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satisfying \l/( 2 )(t , w) = 1. Repeating the steps from (j5.1U|) to (j5.17j) with r = / (r// 2 )(s) (is, 

./to 

provided «?r = zt 2 < 1/e where 

r 2 (t) = jf exp |^ jf B(s') ds'j <fs, (6.16) 

the solution 9^ 2 ^ of the initial value problem ()6.14|) satisfies 

e (2) (M) < — W(-t 2 z) . (6.17) 

T 2 2 

Substituting ()6.4|) and (|6.5|) into ()6.16|) with i = 0, we have 

hm r 2 (0) = /3 hm / 2 _ 3/3/16W = 

to-^-oo t ->—oo L At •W lt " r (sJ 167T — 3/J 

independently of the scale function At, and the power series 

oo 
n=2 

converges if /3 < 167r/3 and \z\ < (er 2 ) ~ . 

6.3 Convergence of the Mayer Series for 47r < j3 < 167r/3 

Q^(0,z) is a candidate to be a majorant of the density /3p(f3,z)/z with the O (z) term omitted 
(see equation (|4.1UJl ). Its linear coefficient C 2 , which satisfies 

C2<^BC 2 + T, t <t<0 (6.18) 

with C 2 (to) = 0, is positive and bounded by 167r/3/(167r — 3/3), in view of the fact r 2 , defined by 
(JHHEJ), solves ()6.18|) as an equality. Hence, by ()3.1|) and ()4.9|) . the Mayer series of Yukawa gas with 
the leading term removed converges uniformly in to for 4tt < (3 < 167r/3 if 

nb n = n lim b n \< C n (6.19) 

A/R 2 ' 

is shown to hold for all n > 2. The inequality (|6.19p . however, does not follow from ()4.3|) because 
the r.h.s. of that equation depends on a 2 (t) = aJ(i)/2 + a2 (t)/2 which, as shown in Subsection 16.21 
diverges as t — > oo. We are going to show that a 2 (t) in the r.h.s. of ()4.3j) can indeed be replaced 
by a>2 (t) and this implies ()6.19|) . 

Define 

- f drU n {r,i) . 

F(t, P n £) = e Js U (s, Pj£; Pj4) ^ c (s, P/£) , (6.20) 
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and suppose n > 3 and \I\ = 2. Using the stability bound ()2.8|) together with ()4.1|) and f)4.2|) . we 

have 



H fdp(^)\F(t,P n ^\ < e^^l^lj] fdp(£, 

lei i£l c l&I 



^2 X i ~ Xj)&j4> C (s, Pl£) 



< e 



n^(s,t) 



r( S )|j c |su P J d 2 x'\ Yl k + ^n^,e,r 



cr'6{-l,l} 



= e " 7(s '* ) r(s) \r\ 4(s) 

by translational invariance. Since a£ satisfies the equation 



(6.21) 



to < t < 



with a£ {t Q ) = (for n = 2 the stability condition (l2~2T| holds with 5 2 = 2), a£(t) < A 2 (t) = C 2 (t)/2 
holds for all t G [to 3 0]- Hence, including the modification (|6.21j) . equation (|4.5j) majorates (|4.3j) 
and (jUTT^|) holds in view of flOJ). 

The Mayer series of Yukawa gas, with the singular part of the norm of ?/> 2 removed, converges if 

16tt - 3/3 



a 16 



and 



z < 



16n/3e 



and this establishes the conjecture up to /3 2 . 



□ 



Remark 6.2 T/ie authors of IBK\I have established convergence of<d for (3 < 167r/3 with O (z) term 
omitted directly from the integral equation L'J.16]) . Our method employs the same ingredients but 
the invariance under translation is used right in the beginning, and not after various manipulations 
of equation \6.21)) . This little detail improves the radius of convergence (see Theorem 4.3 of jBKl) 
and allows us to go beyond the second threshold. 



7 Proof of BGN's Conjecture 

The procedure of Section El will now be extended in order to prove Theorem 12.31 Our proof 
generalizes and includes all results obtained previously. 

The regularized density majorant 0( fc )(/3, z), defined by 0(/3, z) with the singular part removed 
by a Lagrange multipliers, is introduced as follows. If (3 < (3^ = 87r (1 + A;" 1 ) 1 the fc-regularized 
majorant is the solution of equation 

e t = r (^e 2 + z 2 ee z ) + b ( z q z + e-L k ) (7.22) 

where 

L k = L k (t) = 1 + J2 J X ( X - * ( ®z_^Jt, 0) j 

J — 1 \ j— times / 
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with 0(0, z) = 1. We analogously define Q^(t, z) = ty( k '(t,w), where w = f k (t)z and 



k + 1 /"* 
/&(£) = exp <^ — : — / B(s) ds 



t 



(7.23) 



Using 



B~y) w ^w ] + B {^ {k) - 1) - B f 1 - 

* k ' m=2 ^ 



111 — 1 



A 1 



= — k E (™ - * - ijcgv*- 1 < 

m=fc+l 

together with (|7.22j) . is the solution of 



r r 

w 2 ^^„, < —w$> 2 



f; 



fk 



satisfying ty( k \0,w) = 1. Repeating the steps from (j5.10j) to (j5.17j) . the solution 0^" of the initial 
value problem (|7.22jl satisfies 

Q {k \t,z) < —W(-r k z) (7.24) 
provided r k z < 1/e where, by (|6.10p . 



Tk 



t ( k + 1 f f - 
exp < — - — / B(s') ds' ]■ T(s) ds 



(7.25) 



't t ^ Js J P k -/3\K(t) K(t )\K(t)J 

Hence, (|7.24jl defines a convergent power series in z , uniformly in £ , if ft < ftk an d \ z \ < ( er k)~ ■ 
We now examine whether 0( 2r+1 )(O, z) majorizes for ft < ft^r+i the density 

z) = ^ PA - (2b A , 2 z + ■■■ + 2rb A , 2r z 2r ~ l ) , 

z z 

with the leading odd terms up to order 2r — 1 omitted. Because the r.h.s. of equation (J4.3)) depends 
on ci2j(t), j = 1, . . . ,r, which diverges as well as equation (j6.13|) as to — ► oo, inequality (j6.19|) for 
n > 2r do not hold without the modification we have made in the previous section. 

Let F be given by (|6.2(Jj) with n > 2r and |/| <2r and define 



= SU P 



i=2 



«T6{-1,1} B 



0"! 



mt,Pnt)\ 



Note that a® < a n and, by translational invariance, the supreme over x\ is irrelevant. As (|6.21|) 
holds with a 2 substituted by a^, m < 2r, equation ()4.3|) can be replaced by 



a*(*) < j 



1 r* f [n/2] / \ 



m(n-m) a^(s) a n __ m (s) 



n-l 



+ E 

m=[n/2]+l 



// 



m(n -mj a m (sj a;_ m (sj 



(7.26) 
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where [x] indicates the integer part of a real number x. 

We now derive an integral equation for a^'s similar to (|4.3|) . We assume that the improved 
stability condition ()2.21|) holds. Applying the operation 



j „ n n 

— Yldxi^ajvit^x-Xj 

<re{-i,i} n »=i i=i 



I) 



in both sides of ()3.16|) repeating the steps employed for ()4.3|) . with (|2.8|) substituted by (|2.21j) . 
yields 

al(t) <-l da e^^Y(s) J2( n )m(n-m) al(s) a^Js) (7.27) 

1 J m =l V 71 / 

with 5 n > 1/n and n — 1, . . . , 2r. Here, as in equation ([6.21)1 . we have used translational invariance 
of v(t, \x — y\) and, consequently, invariance of t/>£(£) under translation £ — > £' of all components 
= Xj + a by a e R 2 . 

By construction (see analogous equation (|4.4jl ). the Mayer coefficients, regularized so that b 2 j = 
if j <r, satisfy 

„ h < f al/("-l)! if n < 2r 
1 a n /(n-l)! if n > 2r 

and, to conclude the proof, we need to show 

2r ®(+\ 00 (4-\ 

^ (? - 1)! ^ (n - 1) 

j=2 VJ 7 n=2r+l v 7 

which is equivalent to 

r > / aj/(n-l)! = cj if Kn<2r , , 

n - 1 a n /(n-l)! = Cn if n>2r 

where (C n ) n > 2 , given by 

oo 

e w = 1 + c ^ n ' x , 

n 

satisfy, in view of (|7.22jl . the system of equations 



71=2 



C n = -^-(n - l)BC n + — ^ C s C n _j (7.29) 
for n = 2, . . . , k = 2r + 1 and 

C n = nSC n + — C i Cn-j (7.30) 



for n > k, with initial data C n (io) = if n > 1 and Ci(t) = 1 for all t e [to, 0]. 

Comparing (|4.11|) with (|7.22j) . we observe that the Lagrange multiplier L^, k > 1, leads to a 
modest improvement, from n to (k + l)(n — l)/k = n — (A; — n + l)/fc, on the coefficient of the 
linear term C„: 



k — n + 1 1 , , 

n — 1 < n < n . (7.31 

k n 
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On the other hand, equations (|7.27|) and (|7.26J) are equivalent to the system 

„ n-1 

4<(*-^4 + ^£<3<£-i 



i=i 

for n — 2, . . . , 2r + 1 and 

r [«/2] r n-1 

• nl x ^ nl x -» 

Cn ^ TiBc n -\- ^ / C n —j ~\~ ^ / C n —j 

3=1 J=[n/2]+l 

for n > r + 1 with the same initial data (0) = if 1 < n < 2r + 1, c n (0) = if n > 2r + 1 and 
a?(t) = 1, from which, together with (l7~2Hj) . (I?~3UJ) and 5 n < 1/n, equations (l7~2%D follow. 

The Mayer series for the density of the Yukawa gas, with the singular part of the norm of ijj^j 
for j — 2, . . . , r removed, converges if 

P < P2r+i and \z\ < — — 

P2r+lPe 

with f3 n as in (j6.10j) . This concludes the proof of Theorem 12.31 

□ 

We close this section with some comments on the improved stability condition. 

Remark 7.3 The non neutrality constraint is responsible for the improvement \2.21)) since \2.£fy 
saturates if configurations with o~\ + ■ • • + o~ n = are allowed. For n = 2, 

U 2 (t; xi, x 2 , +, +) = U 2 (t; xi, x 2 , -,-)=v (t, \xi - x 2 \) > 

gives 5 2 = 2 and this value exceeds 1/2. For n small, U n (t;£) can be minimized over non neutral 
states £ on a computer. Using the Mathematica software, the minimum for n = 3 

min£7 3 (t;£) > 2min(/i(2w) - 2h(w)) B 

= 2 (h(2w*) - 2h{w*)) B = -2.1162876... ■ B 

is attained for a charge configuration disposed in a straight line with alternate signal separated by 
an equal distance w* = 0.4132466.... This gives a 63 = 0.8837124... larger than 1/3 = 0.3333333.... 

We have also found numerically the minimal configurations for non neutral quadrupoles and 
quintupoles. For n = 4 one of the exceeded charge is expelled to infinity and the remaining charges 
dispose as in the n = 3 case: 

min U±(t; £) > 2 (h{2w*) - 2h(w*)) B = -2.1162876... ■ B 

For n = 5 the minimum is attained for a charge configuration disposed in a straight line with 
alternate signal and symmetric with respect to the plane intercepting the central charge perpendicular 
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to the line: 

minf/ 5 (i;£) > 2 {h{2(w\ + w*)) + h(2w{) + 2h(w\ + w*) 

-2h(2wl + w* 2 ) - 2h(w\) - 2h(w* 2 )) B 
= -4.1879447... • B 

where wl = 0.4399588... and w 2 = 0.2816336... are the distances of the central charge to the 
neighbor and next neighbor charges, respectively. In both case, we have 64 = 1.8837124... > 0.25 
and 5 5 = 0.8 1 20 5 52... > 0.2. 
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